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(54) Clock offset estimation 

(57) A method is provided for estimating the dock 
offset between two clocks located at first and second sta- 
tions (A,B) respectively of a communications network. 
The method involves detenmining for each off a succes- 
sion of update events that provide data concerning the 
offset between the clocks: a first measure e„ which is a 
noisy measure of clock offset; a second measure dn 
indicative of the power of the noise In 0„; and a third 
measure tn indicative of the time interval between the 
update event and a preceding such event. The measures 
derived from each successive update event are used to 



effect a Kalman filtering update process spanning Ksuch 
e/ents. to provide an updated offset estimate. The 
Kalman filtering process may also be used to derive an 
estimate of the error of the clock offset estimate, the 
first time derivative of offset, and an error estimate for 
the latter. A level of noise is preferal3ly injected into the 
Kalman filter to provide the filter with a tracking ability, 
the magnitude of this noise being determined on the 
assumption that convergence has been achieved for the 
offset estimate enror V^. 
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Descriiytion 



Field of the invention 

5 The present invention relates to a method of estimating the clock offset between two clocks located at first and 

second stations respectively of a communications network. 

Background of the Invention 

10 The synchronisation of physically spaced clocks is required in a number of computer and communication applica- 
tions. Thus, for example, If a super-corrputer is to be built by Interconnecting several computers in a network and then 
used to run applications tfiat need a tight clock synchronisation, it is necessary to provide some mechanism for syn- 
chronising the clocks at the individual computers. Another example where clock synchronisation is required is where it 
is desired to analyse a sequence of events taking place across a network where each event has been time-stamped by 

15 a different clock; such a situation might arise in relation to analysing a fault condition on a communication signalling 
network (for example, one operating according to the CCITT Signalling System No. 7 Recommendations). 

Protocols that perform clock distribution across computer networks have been available for a number of years. The 
current lAB standard is the Network Time Protocol, NTP (see, for example, D.L. Mills : "Network Time Protocol (version 
3), Specification. Implementation and Analysis". Network Working Group. RFC 1305, March 1992). The NTP is a general 

20 purpose clock distribution protocol that provides clock synchronisation to arbitrary computer networks and is actively 
used In the Internet. As such, its memory and processing requirements are low. For present purposes, the NTP can be 
considered as made up of a transmit process, a receive process, an update procedure and a local clock process. The 
update process generates the best offset availal>le with reference to a selected reference clock by performing clock 
selection and combination. The local dock part uses the data generated by the update process to adjust the phase and 

25 frequency of the local oscillator with a PLL (Phase Locked Loop) which removes input data noise and provides a stable 
ttmescale between network updates. 

The NTP outputs three quantities: a clock offset with respect to a selected reference clock, an offset error (dispersion) 
with respect to the selected reference clocK and the current round-trip delay to this reference dock. The offset error or 
dispersion can be decomposed in turn into several components: 

30 

1 . - measurement errors due to inherent clock reading enors; 

2. - errors accumulated since the last offset update due to skew enors (here, "skew" refers to the first derivative of 
the offset with respect to time); 

3. - network errors due to delay jitter. 

35 

Normally, delay jitter enors dominate the dispersion budget. The NTP uses a type II PLL to minimise the del^ jitter 
errors. However, PLLs have the problem that even if they work "well" in the sense of good tracking capabilities and 
acceptable errors, it is never known how much "better" one could do. 

It is an object of the present invention to provide an improved dock offset estimation method. 

40 

Summary of the Invention 

According to the present invention, there is provided a method for estimating the clock offset between two docte 
located at first and second stations respectively, the method involving the steps of: 

45 

a) determining for each of a succession of update events that provide data concerning the offset between said 
docks,: 

a first measure 9^ for each update event, which is a noisy measure of said offset; and 
so - a second measure t„ indicative of the relative time interval between the update event and a preceding such event; 

b) using the measures derived from each successive update event to effect a Kalman filtering update process, 
spanning K'such events, to provide an updated offset estimate, the Kalman filtering process being based on the 
following system equations: 

55 

X„,, = X„+f^'„. n^O 
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where X'^ are respectively the dock offset and its first time derivative for the nth update event, and (l4/Jnao 
represents a noise function in the form of a sequence of independent random variables with known upper bound; 
the first time derivative of offset being represented in the Kalman filtering update process as: 

^'n = P n^n 

where: 

Pn = (Pn,7 P^k) 

* (-^n ^n-K+l) 

p^ ^f being determined for each latest update event (for example, by a least squares method) 

IS The system equations are not themselves in a form suitable for application of the Kalman recursions and it is the 
expression of X'n in terms of Pn^n that permits application of the recursions. 

Preferably, the upper bound on the noise is determined for each update event by deriving a third measure d„ in 
respect of that event 

The Kalman filtering process may also be used to derive an estimate of the error of the clock offset estimate, the 
20 first time derivative of offset, and an error estimate for the latter. Advantageously, a level of noise is injected into the 
Kalman filter such as to provide the fitter with a tracking ability, the magnitude of this noise being determined on the 
assumption that convergence has been achieved for the offset estimate error 

In order to remove excessively noisy events from the process, the aforesaid third measure derived from each event 
is preferably compared with a threshold value and if this value is exceeded, the event is discarded. Advantageously, this 
25 threshold value is the average value of the third measure for at least the last preceding update events where M is a 
positive integer. 

The update events used to collect data may take a variety of forms. In one arrangement, each update event is a 
messaging event involving sending an outward message from the first station to the second station with the receipt of 
the outward message by the second station resulting in a corresponding return message being sent from the second 
30 station to the first. In this case, tiie method involves effecting a succession of messaging events, recording for these 
events the timings of message transmission and receipt at the two stations as judged against the corresponding said 
local clocK and determining from these timings for each messaging event: 

said first measure in dependence on the difference between tiie journey times for the outward and return messages 
35 associated with that event; and 

said second measure in dependence on the time interval between the messaging event and the preceding such 
event. 

The messaging events may be effected individually or in groups; in the latter case, for each group, all but the mes- 
40 saging event with the lowest value of the sum of outward and inward message journey times, are discarded. 

Preferably, the time between messaging events (or between groups of such events if the events are sent in groups), 
is set to a minimum value during the accumulation of timings on an initial AC messaging events, and is thereafter deter- 
mined such that the error in offset estimate has grown with time only to a predetermined value before the next messaging 
event or group of such events, as applicable. Is effected, provided that a certain minimum time has elapsed. 
45 Advantageously, for each new offset estimate, a check is made as to whetiier tiie change in estimate between the 
new estimate and the preceding one is within a predetermined limit, the inter-event or inter-group time, as applicable, 
being set to a minimum if this limit is exceeded. 



so 



Brief Description of the Drawings 

A clock-offset estimation method according to the invention will now be described, by way of non-limiting example, 
witti reference to the accompanying diagrammatic drawings, in which: 

. Figure 1 is a diagram showing two stations A and B inter-communicating via a network, each station having its 
55 own clock; 

. Figure 2 is a diagram illustrating the timing of a message event as seen by the clocks at the two stations A and 
B; arKi 

. Figure 3 is a flow chart illustrating tiie control effected over the time interval between tiie sending of bursts of 
messaging events. 
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Best Mode of Carrying Out the Invention 

Rgure 1 illustrates the generalised situation to which the present dock-offset estimation method is applicable. Two 
stations 1 0, herein referred to as stations A and B, inter-communlcate over a network 11 . By way of example, the stations 
5 A and B may be computer workstations and the network 1 1 a LAN (Local Area Network) or a WAN (Wide Area Network). 

Each station A, B has a local clock 12. herein designated clock A for station A and clock B for station B. In order to 
enat>le a correlation of events at the two stations A. B, it is necessary to know the offset of one clock 12 relative to the 
other. This offset will generally change with time and the first time derivative of the offset will hereinafter be referred to 
as the "skew"; the second time derivative of the offset will be referred to as "drift". 
10 The clock-offset estimation method to be described enables an estimation to be made of the offset and skew between 
clocks A and B. Generally, one station (in this case, station A) actively carries out an estimate of the offset of its clock 
A relative to the clock 12 (In this case, clock B) of the other station (station B); station B may be considered as being a 
"time server" and station A as being a client of B. 

This estimation process involves station A effecting a succession of messaging events between itself and station 
IS B. Each messaging event involves the sending of an outward message from station A (source station) to station B 
(destination station) and the return of a return message from station B to station A. Timings associated with at least 
some of these messaging events are recorded and then processed to derive a dock offset estimate. 

In functional terms, station A can be considered as Including, In addition to dock A, five blocks for effecting the 
estimation process. The blocks are a message event block 13 for sending and receiving message-event messages, a 
20 network load monitor 1 4 for discarding "noisy" messaging events, an offset and skew estimation block 1 5, a jump detec- 
tion block 16, and a cooidination block 1 7 for coordinating the operation of the other blocks. Where station A is a computer 
workstation, these blocks will generally be Implemented by appropriate processes run on the workstation processor, tiie 
exception to this being the network interface hardware that conceptually forms part of block 13. 

The functional k)locks 13 to 17 will now be described in turn. 



Messaging Event Block 13 

The messaging event block 1 3 is an-anged to effect a series of messaging events, each involving sending a message 
and receiving the return message generated by station B. The time of sending of the outward message and of receipt 

30 of the return message are noted according to dock A whereas the time of receipt of the outward message and of sending 
back of the return message are noted according to clock B, General mechanisms for providing such a facility are known 
in existing network protocols such as the echoing and time-stamping mechanism of the ICMP protocd associated with 
the TCP/IP protocol suite: with such a mechanism, an outward message is time-stamped on sending by the clock at the 
sending station (station A) and tiie message is then echoed bacK as a return message, from a remote station (station 

35 B) which appends its own time-stamps corresponding to tfie time of receipt and of echo back at the remote station. On 
receipt of the echoed-back message, the original sending station has the required timing information in the received 
message. It will, of course, be appredated that tiie timing Information can be passed from each station A, B to the means 
effecting offset estimation processing (whether or not part of station A) by mechanisms other tiian Including time stamps 
in the outward and return messages. 

40 In the present exanrple. it is assumed that the stations A. B run the ICMP protocol and use its mechanisms to effect 
and time the messaging events, station B providing timestamps for botti message receipt and echoback. 

Referring now to Figure 2, the flow of a messaging event (taken as the n*^ such event) is illustrated relative to the 
time lines 20, 21 for clocks A and B respectively. Thus, an outward message (represented by line 22) is first sent from 
station A to station B which then sends back a return message (line 23) to station A. The time of sending of the outward 

45 message 22 is recorded as being at time r„ by clock A (the suffix "n" indicating that tiie timing relates to an n^*^ messaging 
event); note that in practice, time r„ is sllghtiy k>ef6re the outward message is actually transmitted from station A, this 
interval being a handling interval needed by station A to effect ti-ansmission after timestamping at 7^. The receipt of tiie 
outward message 22 at station B is recorded by a time stanrp at time Rp according to dock B (note that for clarity, timings 
according to clock B are primed); between receipt of the outward message 17 and /?„ is a small handling interval. The 

so sending by station B of the return message 23 is recorded by a time stamp at time S'„; between S'„ and transmission 
of the return message 23 is anotiier small handling interval. The time of receipt of the return message at station A is 
recorded as time according to dock A, there again being a small handling interval between receipt of the return 
message and time A^ The messaging event block 13 thljtf has available time values 7„, S* „. A„ following the n^ 
messaging evem (note that If stations only provides one time starnp for message receipt ardechobac^^ R'n = 

ss Generally, of course, tiiere will be an offset between clocks A and B. Thus, when dock B records times R'n* S'n. 
clock A will show different values respectively. More particularly, if X(t) is tiie dock offset at source time t (clock 
A time), then: 



25 



R„^R„^X(RJ 



(1) 
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s„=s'„+xrsj (2) 

Similarly, times >A„ and 7";, at the source (station A) correspond to times A'^ and T„ at the destination (station B). 
For the n^ messaging event, block 13 calculates the quantities: 

Half round trip delay c/„ = g (3; 

10 PacketOffset e„ - ' ^'"^g ' ^"^ (4) 

The quantities d„ and e„ are be directly computed from the timestamps. 

The larger d^, the "noisier" the offset information of the n*^ messaging event. In order to reduce the noise (that is. 
75 to reduce the half round trip delay og, rather than sending Isolated messaging events, block 13 is anranged to send 

bursts each of several messaging events, with, for instance, the events in a burst being separated by about one second. 

The first messaging-event packet in a burst acts as a "warning shot" and has the effect of moving the destination address 

into the cache of routers in the network. As a result, subsequent messaging-event packets of the same burst normally 

end up with a smaller half round trip delay. For every burst of messaging events, block 13 selects the messaging event 
20 with the smallest value of d„ 

Thus, for the n*^ burst of messaging events, the block 13 selects the messaging event with the smallest value of 

and treats this selected messaging event as the n^ messaging event; block 13 then proceeds to output the following 

values for this selected event: 

Network l^ad Monitor 14 

The network load monitor 14 tracks the network load. The network load is estimated from the half round trip delay 
30 values. Upon output of the n*^ messaging event data by block 1 3. the load monitor 1 4 first uses the half round trip delay 
value dn to update an admission threshold A to a new value Preferably, A„ is computed as the average half round 
trip delay seen over the last M selected messaging events from the last M messaging-event bursts. It is the value of M 
rather than which average (e.g. arithmetic mean, median,...) that plays an important in monitoring the network load, a 
large value of M making the monitor less adaptable and a small one making the monitor redundant {M = 1 would mean 
35 that each messaging event would always be accepted). 

After the threshold value A^ has been updated, the value of cf„ is compared to the threshold. If d^ is less than A„ 
then the n**^ selected messaging event is considered valid and the values output from block 13 are passed to the esti- 
mation block 1 5 for dock offset and skew estimation; otherwise, the n^^ messaging event is discarded and the coordination 
block 1 7 informed. In this latter case, the coordination block 1 7 initiates a new messaging event Ixjrsl by the block 1 3. 



40 



45 



55 



Offset and Skew Estimation Block 15 

Block 1 5 estimates the dock offset X(t) and the clock skew X'(t) from the valid selected messaging event observations 
{dfi. Tf^ Af^, nsO. This is carried out using Kalman filtering as described below. 



Clod< Offset Model Assumptions and Notation 
Evolution Equation: 

50 It is assumed that the clock offset X(t) at source time t has an evolution equation of the form: 



X(t)^x-i^jX'(s)ds (5) 



the quantity X(0) being called the initial dock offset and the function X'() being the dock skew (first time derivative of 
offset). 
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Assumptions: 

(1) - the drift (that is, the second time derivative of offset) is practically zero so that between two consecutive valid 
messaging events, the clock skew remains practically constant; 
5 (2) - the clock skew per unit time is assumed to be very small (a few parts per million): 

\X'(t)\ «1. ^sO 

(3) - protocol stack processing and routing delays (source-destination: destination-source) are on average symmet- 
10 rical. 

Assumption (3) is often not met in practice; however, this is not critical as this only leads to a constantly biased offset 
value, the bias being dependent on the degree of asymmetry. In contrast, assumptions (1) and (2) are more critical in 
nature (but see later for how a failure of assumption (1) can be handled). 

15 

Notation: 

Let: 

Basic Equations 

25 

Singie-Event Equation - The first equation that can be set down for relating messaging event data to the clock 
offset nrKxlel. concerns each event taken in isolation. The quantity (which is half of the difference between outward 
and inward packet journey times) would actually be a true measure of clock offset were it not for the variations from the 
symmetry assumptions noted above that affect any single messaging event. It may therefore be said that: 

en = X„+ M^„,naO (6) 

where (W/^Jnao represents a sequence of independent variables of zero mean (if assumption 3 holds) bounded by 
and witii variance {a^^, W„ is called the measurement noise. The initial value Xo \sa random variable of mean Oq and 
35 variance (ao)^. 

Event-Sequence Equation - The second equation that can be derived for linking the measurement results with tiie 
offset model, concerns changes taking place from one messaging event to the next. Clearly, tiie offset at event (n+1) 
will be related to the offset at event (n) by the dock skew (first time derivative of offset); however, there will also be some 
randomising factor involved. The second equation has the form: 

40 

X„^, = X„ + f„X„+ W'„^,.n^O (7) 

where: 

and (^Vn)n£i represents a sequence of independent zero mean random variables with variance (v^. IV'„ Is called the 
driving noise. The random variable sequences and are independent of each other. 

so Kalman FiHerinq 

Based on equations (6) and (7) above, a Kalman filter is provided to recursively confute the clock offset estimate: 

^ ^n = ElXn|eo,.-.,0n], « > 0, 
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and the mean square error: 



In this process, the clock skew (the first time derivative of offset) is estimated from the measurement data. More generally, 
there may be situations where it is important to capture the offset dependency on the skew, the drift and higher derivatives 
10 (this would permit assumption (1) above to be breached); in the following it is shown how this can be done with a K 
dimensional Kalman filter though, for darity. It is assumed that only the skew is needed (the same approach being usable 
for higher order time derivatives). 



Offset Estimate - First, the skew is written as: 



IS 



20 



where: 



Pn'=(Pn,r Pn,K) 



(9) 

(10) 
(11) 



The system equations (6) and (7) can now be rewritten in the following form (superscript Tdenoting the transposition 
ss operator): 



30 



where: 



H=(1,0 0)^ 



(12) 
(13) 

(14) 



35 



40 



45 



/ 1 + a„,i q;„,2 
1 0 



JPn = 



0 

0 
0 



0 
0 



0 0 



0 
0 

1 



0 

0 
0 



(15) 



50 



M^'^i=0^n.n.O 0)^ 



(16) 
(17) 



denotes the estimate of Xn given e o e^,. and the matrix Cn denotes the Kx Kcovariance matrix of Xn. that is: 

55 then the theory of Kalman filtering provides formulas for the updates of Cn and Xn at time An^^p 
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These updates are as follows: 



10 
15 

where: 

20 



(19) 



C -(FCf'^Q ) (^nC„F„"^-^Qn.,)HHV„C„F/^Q^,) 



= H^C„H (21) 
«„' = a„C„a/ (22) 



25 O„^,=diag{t/„^0 0} (24) 

The initial values are set after the first K messaging events as: 

A 



30 



ATo — (6^+i>-»->9a)^ (25) 



35 Co = diagUa.K^^,)'' (a^)^} (26) 

Qq = cliag{(i^o)^.0 0} (27) 

The offset estimate at time fA„^y + 7"„^,;/2 is: 

^ ^ ... p ^ ^ 

•„ + (^„*,) (e ^ , - (X „ + f „X- „)) n a 0 

and its mean square error is: 

45 (V„^,)*-(a„^,)V„*,)^/)sO (29) 

where: 



Skew Estimation - To enable the dock offset to be derived according to the above formulae, the vector p„ is estimated 
55 using a least squares estimate. If E = (e ^) is the {K- 1)xK matrix given by : 

1if/ = y (31) 
e^ = -1 if y= /+7 
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0 otherwise 

where 1^^{K-r) and ^sMK, then the (K-l) vector with i-th component, 

n-i+l ■ ^ n-i)^n-i» 

1 s^K^can be written as T„X„ Where T„ = diag{1/fn.7, .... l/^n/c+tl^n- 
The covariance matrix G„ of the random vector T„X„ is equal to: 

G„=7-„C„r„^ (32) 
The function to minimise in order to obtain a least squares estimate of p and of is: 

= (TSn-X\efGjiTX'X\e) (33) 



Where is the generalised inverse matrix of G„ and e = (1 1)^ From this it can be shown that the clock skew 

20 estimate is then: 



A / \A 



(34) 



and thus. 



p ^l!^nlln (35) 
e G„ e 



and the skew variance is given by: 



(a^,)^=Var(Xg« ' ^ (36) 



Dynamic Behaviour - There is a trade-off between the Kalman filter tracking ability and the Kalman filter en-or 
estimate. The driving noise power (v>^ is used to control this trade-off. If the driving noise is too large. V^^ tends to a 
value that is unneeessarily large. The question to address is what is the right power of driving noise that needs to be 
45 injected into the Kalman filter to preserve its tracking ability without creating error estimates that are unnecessarily large. 

The Injection noise le^el used is that which would correspond to having achieved enror convergence, that is, 
V^^ V^n+ 1 • "^^^ at reception of the (n+1)th messaging event, the level of noise to be injected into the Kalman 
filter is given by: 



so 



ss where [x]* = max {0.x} . 

This approach seems to enhance the tracking capability of the filter at the expense of the offset en-or performance. 
This is justified since the errors obtained by the Kalman f flter are normally an order of magnitude below the half-round 
trip delay. 
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Implementation 



For the (n+1)th valid messaging event, the estimation block 15 has the following new data to work with: 

^ . the half round trip time which provides an upper bound on the measurement noise; it is generally used to 
5 determine the standard deviation of the measurement noise when evaluating the above formulae: 
1 - the packet offset; 
tn ' the inter-event time (derived according to equation 8 above). 
This data is then used to update the offset estimate as follows: 

10 1 ) the matrix F„ is determined by: 

calculating the vector p„ using formula (35); 
- calculating / using formula (17); 
calculating F„ using formula (1 5); 

2) the matrix ^ is determined by: 



15 



calculating the driving noise power ( 7)2 using formula (37), the value of V„ having been determined in relation 
to the previous messaging event; 
20 - calculating / using formula (24) ; 

3) the value of is determined using formula (22), the value of C„ having been determined in relation to the 
previous messaging event; 

4) the offset is estimated using formula (19); 

25 5) the mean square error of the offset is estimated by: 

determining the value of ( y)2 using formula (30); 
determining the mean square error using formula (29); 

30 6) a skew estimate is derived using formula (9) together with an estimate of its mean square en-or using formula (36) ; 
7) the value of C^.^ is determined using formula (20). 

It will be appreciated that the numerical computations of the /^dimensional Kalman filter need to be done carefully 
in order to avoid numerical problems (such as division ty zero); appropriate techniques for avoiding such problems will 
35 be apparent to persons skilled in the art. 

Jump Test Block 16 

A highly unexpected value of e„ may indicate a drastic change in the network behaviour or conditions; in particular. 
40 it may indicate a jump in the value of clock offset. Clock offset jumps will generally be rare and much larger than fluctu- 
ations due to noisy data. Block 16 tests for such Jumps and if one Is detected, block 16 inlbrms the coordination block 
16 which sets the inter-burst time to minimum to ot>tain further data as quickly as possit)le. 

The jump test carried out by block 16 also has the function of ensuring that the Kalman filter maintains its tracking 
akdilities - if it does not. then this will be picked up by the junrp test. 
45 The nature of the jump test is as follows (this being one example only, other tests being possible). On the k>asis that 
the clock skew will only vary slowly with time, the hypothesis: 

^0''^n^1-^n+t„X'„ (38) 

SO against the hypothesis: 

If hypothesis IHo holds, then it is taken that no offset jump has oocun^ed. Under hypothesis Ho. it may be expected that: 

e„^,- H^(F^„+ W'„,,)^W„,, (40) 
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and since 



n./ = e,.,; - ti'FjL, (41) 



TO and Xr, are independent, it follows that the quantity: has zero mean and unit variance; hypothecs Ho can be rejected if, 

for example, I ^ ^ I > z , where z is a pre-specified value. 

The numerator in formula (41) the absolute size of the clock offset jump. Thus if one is interested in detecting 

jumps of a fixed size, regardless of the normalisation introduced by the denominator of formula (41 ). It suffices to compare 

the numerator in formula (41) with the fixed size threshold off interest. 
15 If hypothesis Ho is accepted no action is taken. However, if hypothesis Hq is rejected, the coordination block 1 7 is 

informed. 



CiQQrdination Block 17 

20 The coordination block 1 7 coordinates the operation of the other blocks and, in particular, controls the time between 
message-event bursts in dependence on how well the estimation process is doing according to feecfl>ack provided by 
the other blocks. Clearly, at start up when data on Kvalid selected messaging events is required to start off the Kalman 
filter, the inter-burst time is preferably set at a minimum; however, once a steady state condition has been achieved, 
then' the inter-burst time wfll depend on the allowable level of en'or for the offset estimate (the longer the time between 

25 bursts, the greater the potential error). More particularly, the offset error f time units after time (4 „ + 7^)12 is given by: 

V^'CO = (H + /P ^C„(H ^ + J ^ (42) 

where: 

30 

^^(0)= V„^ = h'^C„H (43) 

is the offset error at time ( A „ + T„)/2 

The offset error increase will typically be controlled not to exceed a value ec/„ where e is a preassigned quantty. If 
35 {e^d„^' V„^)>0, the (n + 1)th burst will be sent Mime units after the nth burst, where fis the positive solution of: 

(a^.)V + 2W.(E^d„^-V/) (44) 



Since usually 

and since also there will be a minimum allowable time between bursts, the inter-burst time f will be set such that: 
- if (e^d„ ^ - V„^) > 0 . f is the greater of: 

(^n g • ) 



where 

40 



55 

and 

(mwi 

- othenivise. t = f ^, 
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Figure 3 is a flow chart illustrating the operation of the coordination block 1 7 in controlling the sending of message-event 
bursts by block 1 3. At startup, a counter has its value c set to zero (step 30); this counter is used to count the number 
of accumulated valid messaging events (Ksuch events being required to start off the Kalman filter). Next, the inter-bust 
time / is set to fmin (st^P ^^^^ ^i^® ^ 3^) ^® coordination block 1 7 triggers the block 13 (step 33) to send 

5 a messaging event burst. 

If block 13 receives no return message from station B then this Is reported to the coordination block 1 7 and block 
1 7 returns to step 31 (see step 34). 

If. however, at least one return message is received by block 13 for the messaging-event burst, block 17 is informed 
that a messaging event has been effected and block 1 7 then proceeds to check the count value c (step 35) If c is less 
10 than the required number K needed to run the Kalman filter, tiie value of c is incremented (step 36) and block 1 7 returns 
to step 32; however, if c equals K, then block 1 7 awaits a "valid event" indication from the network load monitor 1 4. 

If the network load monitor 14 reports that no valid messaging event results from a burst (tested in step 37) , then 
block 17 returns to step 31. However, if a valid messaging event is reported, block 17 next awaits receipt of the new 
offset skew and error estimates from the estimation block 15 (step 38). 
15 Thereafter, if block 16 reports that no offset jump has been detected (tested in step 39), block 15 proceeds to output 
the offset estimate (step 40) and determine a new value for / (step 41) before returning to step 32. In the presence of 
an offset jump, block 1 5 returns to step 31 . 

It should be noted tiiat the offset estimate determined by block 15 is the offset at time (A^-^ r„ )/2 . The offset 
estimate at time is determined by taking account of the offset change due to the value of the skew over the period 
20 {A„- T„)/2 . The offset at any intermediate time between two estimates can be similarly derived. 

Variance Estimation 

In the foregoing, the variance a„2 of the measurement noise was equated to its known upper bound dn^. However. 
25 the quantity VVt introduced in formula (41) must have a variance of 1 which is not achieved with a„^,sd„^, because 
i is too high, it is, in feet, possible to define a parameter y. between 1 and 0 in value, such tiiat if a ^ , - yd„ ^ ^ ttnen 
the variance of Y^^ is equal to 1. Furthermore, in this case, all variances are scaled by y: as a result, the offset en^or 
\/„2 which was previously derived using a„ ^ =d„^. now becomes scaled down to V2\/„2. 

The value of y is estimated for the ATth valid messaging event as tiie running variance of all the quantities Y„ to date, 
30 where tiie quantities ( comprise the sequence of values obtained from formula (41) witti d^. 

Vgriarrtg 

It will be appreciated ttiat many variants are possible on the foregoing. For example, the steps taken to exclude 

35 noisy messaging events may be omitted, though again this is not preferred. 

Although in the described embodiment, messaging events liave been used to collect the measurement data to which 
the Kalman filtering process is applied, it is possible to collect tiie necessary update data in other ways. For example, 
one station might periodically inform the otiier of what it thinks its offset is; alternatively, in appropriate circumstances 
one could infer the offset of a station by the way it generated traffic. The update measures required are an offset measure 

40 e„ (which will inherentiy be noisy), an upper bound dn on the noise of the offset measure, and a knowledge of the time 
between successive update measures. With regard to tiie upper bound on tiie noise power in the offset measure, if tiie 
environment concerned is sufficientiy well known, it is feasible to dispense with making an upper bound measure for 
each update and simply use a predetermined value; however, this is not preferred. 

45 Claims 

1. A method of estimating the dock offset between two clocks located at first and second stations respectively of a 
communications network, said method comprising the steps of: 

so a) determining for each of a succession of update events that provide data concerning the offset between said 

clocks,: 

a first measure 6^' for each update event, which is a noisy measure of said offset; and 
a second measure tn indicative of tiie relative time interval between tiie update event and a preceding such 
55 event: 



b) using the measures derived from each successive update event to effect a Kalman filtering update process, 
spanning K such events, to provide an updated offset estimate, the Kalman filtering process being based on 
the following system equations: 
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where X„. X'„ are respectively the clock offset and Its first time derivative for the rth update event, and ( M/^neO 
represents a noise function in the form of a sequence of independent random variables with known upper bound: 
the first time derivative of offset being represented in the Kalman filtering update process as: 

where: 

Pn = (Pn,7 Pn,/c) 

IS ^n-^^n 

p„ ;f being determined for each latest update event. 

2. A method according to claim 1, wherein said upper bound on the noise is determined for each update event by 
20 deriving a third measure </„ in respect of that event. 

3. A method according to claim 2, wherein said Kalman filtering process is also used to derive an estimate of the en-or 
V„ of the dock offset estimate. 

55 4. A method according to claim 3. wherein a level of noise is injected into the Kalman filter such as to provide the filter 
with a tracking ability, the magnitude of this noise being determined on the assumption that convergence has been 
achieved for the offset estimate error V„ 

5. A method according to daim 2. wherein an estimate is also made of the value off the first time derivative of the offset 
30 by using covariance data obtained from the Kalman filtering process. 

6. A method according to daim 1 . wherein between offset estimates, the actual current dock offset is estimated by 
combining the latest offset estimate with a component dependent on the lapsed time since that estimate and the 
first time derivative of the offset. 

35 

7. A method according to claim 2. wherein a said update event is discarded if the said third measure derived therefrom 
exceeds a predetermined value. 

8. A method according to daim 7. wherein said predetermined value is the average value of said third measure for at 
40 least the last M preceding update events where /W is a positive integer. 

9. A method according to any one of the preceding claims, wherein each said update event is a messaging event 
involving sending an outward message from said first station to said second station with the receipt of said outward 
message by the second station resulting in a oon-esponding return message being sent from the second station to 

45 the first; the method involving effecting a succession off said messaging events, recording for these events the timings 
of message transmission and receipt at the two stations as judged against the con-esponding said local dock and 
determining from said timings for each messaging event: 

. eaid first measure in dependence on the difference between the journey times for the outward and return mes- 
so sages assodated with that event; and 

. said second measure in dependence on the time interval between the messaging event and the preceding such 

event. 

1 0. A method according to daim 9. wherein said messaging events are effected in groups each of at least two messaging 
55 events, and for each group, all but the messaging event with the lowest value of the sum of outward and inward 

message journey times, are discarded. 

1 1 . A method according to claim 9 or daim 1 0. wherein the time between messaging events or between groups of such 
events if the events are sent in groups. Is set to a minimum value during the accumulation of timings on an initial K 
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messaging events, and is thereafter determined 
such that the en^or in offset estinnate has grown with 
time only to a predetermined value before the next 
messaging event or group of such events, as appli- 
cable, is effected, provided that a certain minimum 
time has elapsed. 

12. A method according to claim 11. wherein for each 
new offset estimate, a check is made as to whether 
the change in estimate between the new estimate 
and the preceding one is within a predetermined 
limit, the inter-event or inter-group time, as applica- 
ble, being set to a minimum if said limit is exceeded. 
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